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NOMENCLATURE 


a wing  lift  curve  slope  3L/ 3a 

A ratio  of  rms  incremental  load  factor  response  to 

rms  vertical  gust  velocity 

AR  wing  aspect  ratio 

c reference  chord  of  wing  or  vertical  tail  depending 

on  subscript 

e natural  logarithm  base 

f unsteady  lift  function 

f.  component  of  unsteady  lift  function  due  to  aircraft 

motion  defined  in  the  frequency  domain 

f2  component  of  unsteady  lift  function  due  to  gust 

1 penetration  defined  in  the  frequency  domain 

F,G  Theodorsen  functions  for  unsteady  lift 

g gravitational  acceleration 

h altitude 

H transfer  function  (see  Eq.  (6)) 

k nondimens ional  temporal  frequency,  wc/2U 

k upper  limit  of  integration  for  the  calculation  of 

c A and  NQ 

k exceedance  function  from  Ref.  1 

o 

In  natural  base  logarithm 

L characteristic  scale  of  turbulence 

M Mach  number 

N number  of  times  a given  load  level  is  exceeded 

P ratio  of  flight  time  in  heavy  turbulence  to  total 

flight  time 

S wing  reference  area 

t time 
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V 


T 


'I 


u 

w 

W 

x 

X 

z 


mean  speed  of  aircraft  relative  to  air 

vertical  component  of  fluctuating  velocity  of 
aircraft  relative  to  air 

weight  of  aircraft 

position  coordinate  parallel  to  U and  positive 
in  the  forward  direction 

integration  limit  for  x position 

position  coordinate  in  vertical  direction  positive 
upwards 


a 

An 

y 

p 

o 

k 


* 


w 

fl 


angle  of  attack 

incremental  load  factor  ~ (L  - W)/g 
aircraft  density  parameter  2W/apcgS 
density  of  air 

standard  deviation  or  rms  value  of  a fluctuating  function 

power  spectral  density  function 

gust  alleviation  factor  from  Ref . 1 

temporal  frequency 

spatial  frequency  of  turbulence 


Superscript 

(*)  differentiation  with  respect  to  time 

Subscripts 

w vertical  fluctuating  velocity  component  or  wing 

An  incremental  load  factor  due  to  a vertical  gust 

VT  vertical  tail 

d desired 


vi 
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1 . INTRODUCTION 


This  report  describes  a simple  technique  for  estimating  loads 
on  aerodynamics  surfaces  produced  by  atmospheric  turbulence.  It 
is  intended  to  provide  the  designer  with  a technique  for  estimating 
gust  loads  more  realistically  than  the  discrete  gust  method  and 
with  little  increase  in  complexity. 

At  present, aircraft  covered  by  FAR23  are  designed  to  safely 
encounter  a discrete  gust  with  specific  size  and  shape.  However, 
the  continuous  and  random  nature  of  atmospheric  turbulence 
suggests  that  the  designer  consider  not  only  the  magnitude 
of  the  load  but  also  its  frequency  of  occurrence.  A typical 
histogram  of  load  exceedances  such  as  shown  in  Fig.  1 would  allow 
for  a more  rational  choice  of  gust  design  load  and  also  provide 
data  relevant  to  fatigue  studies.  For  example,  aircraft  struc- 
tures could  be  designed  to  handle  a gust  load  level  which  is  ex- 
pected to  be  exceed  only  10  times  in  30,000  hours  of  flight. 

Using  Fig.  1,  this  would  correspond  to  a 10,000  pound  load. 

There  are  several  methods  which  can  be  used  to  predict  the 
load  level  exceedance  distribution  for  a particular  aircraft 
and  series  of  mission  profiles.  Estimates  can  be  obtained  using 
either  time  or  frequency  domain  (power  spectral)  representations 
of  atmospheric  turbulence  and  aircraft  response.  The  complexity 
of  these  methods  and  their  ease  of  application  depend  mainly  on 
the  complexity  of  the  aircraft  gust  response  model.  References 
1-7  develop,  describe  and  evaluate  several  load  exceedance  cal- 
culation techniques  in  both  the  time  and  frequency  domain. 

This  design  guide  illustrates  a method  which  uses  the  power 
spectral  representation  of  atmospheric  turbulence  and  a one  degree 
of  freedom  model  of  aircraft  gust  response.  It  is  easily  ac- 
complished using  several  design  charts  and  is  therefore  well  within 
the  capabilities  of  all  aircraft  designers. 


Number  of  times  a given  magnitude  of  incremental 
wing  load  is  exceeded  in  30,000  hours 


The  methods  outlined  here  should  be  restricted  to  use  with 
fairly  rigid  aircraft  flying  well  below  transonic  flight  speeds. 
Also,  the  small  perturbation  lateral  and  longitudinal  responses 
should  be  well  damped  and  should  not  be  strongly  coupled. 

Section  2 details  the  design  technique  and  gives  the 
mathematical  development  of  the  design  charts.  Section  3 gives 
a step-by-step  description  of  the  design  procedure  along  with 
an  illustration  for  a specific  aircraft  and  flight  condition. 

For  a more  detailed  comparative  evaluation  of  this  design  proce- 
dure and  other  power  spectral  techniques  the  reader  is  referred 
to  Ref.  7. 
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This  section  is  intended  to  provide  the  designer  with  a 
basic  understanding  of  the  one-degree-of-freedom  power  spectral 
gust  load  design  procedure.  It  briefly  reviews  the  power  spectral 
approach  to  the  problem,  modeling  the  atmospheric  turbulence, 
modeling  the  aircraft  gust  response,  and  finally  determining  load 
level  exceedance  curves.  No  attempt  is  made  to  compare  this 
technique  with  any  other  method  . The  technique  is  developed 
for  specific  application  to  gust  loading  of  the  wing  but  can 
also  be  used  for  other  aerodynamic  surfaces  as  will  be  pointed 
out  later. 

Assume  that  an  aircraft  flying  straight  and  level  with  air- 
speed U encounters  a patch  of  turbulence  characterized  by  a 
vertical  velocity  w(x)  (Fig.  2).  Th j forward  velocity  of 
the  aircraft  will  make  traversing  the  turbulence  patch  appear 
like  a time-varying  velocity  w(t)  . If  the  vertical  velocity 
is  upwards  it  will  increase  the  lift  on  the  aircraft  wing  and 
will  cause  the  aircraft  to  accelerate  upwards,  z(t)  . The  accel- 
eration z is  usually  expressed  as  the  incremental  load  factor 
An  = z/g  . It  is  this  load  factor  response  which  must  be 
estimated  and  used  in  the  structural  design  specifications  of 
the  aircraft. 

Due  to  the  random  nature  of  atmospheric  turbulence,  it  is 
easier  to  work  with  the  spatial  frequency  distribution  of 
vertical  velocity  (Power  Spectrum)  4>w(^)  , rather  than  w(x) 
or  w(t)  These  functions  are  related  by 

X/2  2 

f w(x)eifi  dx  (1) 

-X/2 

where  n is  the  spatial  frequency.  For  this  and  most  gust  load 
estimation  techniques,  the  Von  Karman  power  spectral  representa- 
tion of  atmospheric  turbulence  is  used 


W X—  X 


A 
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(2) 


1 + | (1.339Ln)2 

*w<U2)  it  £i  + (i.339LJ2)2]11/6 


where  a is  the  root  mean  square  turbulent  velocity  and  L is 
w 

the  scale  of  turbulence.  Reference  1 recommends  that  a value  of 
L = 750  ft  be  used  for  design  purposes.  Reference  9 recommends 
ow(h)  for  severe  turbulence  resulting  from  strong  convection  or 
storms  as  shown  in  Fig.  3.  This  spatial  frequency  distribution 
appears  to  an  aircraft  traversing  it  as  a temporal  frequency  distri- 


bution. 


*w<u)  " u »w<L“> 


(3) 


in  which  to  = Uft 

The  basic  physical  model  of  aircraft  response  to  turbulence 
is  that  an  upward  gust  velocity,  w , produces  an  acceleration,  z 
in  the  vertical  direction. 


l z(c)  - 


apU2S 


f (w,t) 


(4) 


where 


W = aircraft  weight 
g = gravitational  constant 
a = wing  lift  curve  slope 
p = density  of  air 
U = mean  airspeed 
S = wing  reference  area 
f (w,t)=  unsteady  lift  function 

Expressing  the  normal  acceleration  as  incremental  load  factor  An 

An(t)  = f(w,t)  (5) 

This  equation  is  transformed  into  the  frequency  domain  and  expressed 
in  transfer  function  form 
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•*» 


*An<k>  « ^ 


f1(k)f2(k) 


"1 1/2 


(6) 


where  the  unsteady  lift  function  is  divided  into  the  aircraft 
motion  effect,  f^  , and  the  gust  penetration  effect,  f2  • 


where 


fx(k) 


16u2 


f2(k) 


(¥)  ♦ (*v * m 


1 - M2 

1 - M2  + 1.5wk  + ir2Mk2 


(7) 


(8) 


2W 

apcgS 


(aircraf  : density  parameter) 


k = ^ (nondimens i< nal  frequency) 

and  F , G are  the  Theodorsen  functions  used  with  two-dimensional 
airfoils  in  oscillatory  incompressible  flow.  For  design  purposes 
f2  is  evaluated  at  M = .2  which  should  be  satisfactory  for 
0 < M < .5  . 

The  output  power  spectra  (k)  is  computed  using 


Wk> 


Hin(k)  2 *„(k) 


(*#0  fi<k»2<k>*w<k> 


The  r.m.s.  value  of  the  incremental  load  factor  is 

00 

11/ 2 


°An  = J<t»An(k)dk  * TJT  Vw 


where  from  Eq.  (9)  and  Eq.  (10) 


(9) 


(10) 


r 


♦„00 


f1(k)f2(k)  -^7 


dk 


w 


1/2 


(ID 


Note  that  the  upper  limit  of  integration  has  been  changed  from 

TT 


infinity  to  kc  which  for  numerical  purpose  is  set  to  kc 


The  structural  response  parameter,  A 


of 


An 


An 


to 


w 


is  therefore 
o 

’w 


A _ °An  _ apsu  v _ U v 
AAn  " o„  TUT  % peg  % 


which  is  the  ratio 


(12) 


K 


4> 


is  a function  of  y (from  f^)  and  2L/c  (from  <j>w)  and  is 
shown  in  Fig.  4.  Using  the  results  contained  in  this  figure  plus 
Eq.  (12)  the  response  parameter  A.  can  be  evaluated. 


An 

The  rate  of  zero  crossings  for  An 


is  given  by  the  relation 


„ _ 1 An 

N° ' 


kc 

lU  ( 

2* 


"in  *w<k)dk 


in  the  positive  direction 
]l/2 


1] 


(13) 


An 


This  can  be  written  as  a function,  kQ  , which  is  given  by  y and 
2L/c  (Fig.  5)  and  a velocity  term 


N - JL  k 

O TTC  O 


(14) 


Using  Eq.  (13)  and  Fig.  5 the  zero  crossing  parameter  Nq  can  be 
evaluated.  Finally,  for  severe  loads,  the  number  of  times  per  second 
that  any  given  loading  An^  is  exceeded  is 


N(And)  = PNQe 


- An,/A.  a , 
d'  An  w 


(15) 


P is  the  ratio  of  flight  time  in  severe  turbulence  to  total  flight 

-4 


where 

time.  Reference  9 recommends  P = 5 v 10-M  at  altitudes  below 
40000  ft.  For  altitudes  above  40000  ft  consult  Fig.  6 for  the 
recommended  value  of  P . 
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Gust  alleviation  factor 


In  addition  to  wing  loads  it  is  also  desirable  to  determine 
design  values  of  the  gust  loads  on  the  horizontal  and  vertical  tail 
surfaces.  In  Ref.  7 it  was  found  that  a zero  degree  of  freedom 
(ODF)  model  produced  consistent  and  useful  estimates  of  gust  loads 
in  the  vertical  tail.  Estimation  of  gust  loads  on  the  horizontal 
tail  requires  the  use  of  the  more  complex  two-degree-of-freedom 
modeling  since  it  is  significantly  influenced  by  the  downwash  from 
the  wing.  Alternatively,  use  the  discrete  gust  design  procedure. 

The  ODF  model  is  a subcase  of  the  IDF  model  already  discussed. 
It  is  equivalent  to  the  IDF  model  with  infinite  aircraft  weight. 

In  this  case,  however,  the  designer  is  interested  in  vertical  tail 
load  expressed  in  units  of  force  (lbs)  instead  of  units  of  accelera- 
tion. In  this  case  the  r.m.s.  value  of  vertical  tail  load,  , 

is  given  by  the  relations 

a — apSU  „■ 
avt  ~~l  % 


VT 

(16) 


where 

\ 

* 


aPyT  " ^VT°w 


(17) 


a = vertical  tail  lift  curve  slope 
S = vertical  tail  area 


is  obtained  from  Fig.  4 using  y = «>  and  the  mean  aerodynamic 
chord  of  the  vertical  tail  to  calculate  2L/cv^  = 1500/c^  . 

The  evaluation  of  N proceeds  by  determining  k at 

°VT  ° 

y = 00  and  2L/c  = 150Q/cVT  from  Fig.  5 and  substituting  into  the 
relation 


(18) 


As  shown  in  Eq.  (15),  the  estimation  of  A and  NQ  for  a 
particular  aircraft  and  flight  conditions  completely  specifies  the 
rate  at  which  any  given  load  level  is  likely  to  be  exceeded 
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* 


' 


However,  if  the  design  criterion  specifies  the  rate  of  exceedance, 
N , and  the  corresponding  gust  design  load  factor  is  to  be  deter- 
mined, a more  useful  form  of  this  relation  is 


/ n \”^ 

And  = AAn°w  £n 


(16) 


The  equivalent  form  of  this  equation  for  use  in  specifying 
a vertical  tail  design  load  is 

v-1 


PVr  AVT°y/j 


Jtn 


N 


PN 


(17) 


VT 


It  is  also  possible  to  obtain  design  loads  for  the  wing  and 
tail  in  a simpler  fashion  analogous  to  the  discrete  gust  design  method. 

K 


ud 


And  " AAn  Kg 


Where  = gust  design  velocity  (usually  50  fps . ) this 


o , or  P,  but  it 
w’  ’ 


method  does  not  require  an  evaluation  of  Nq , 
also  does  not  yield  any  frequency  of  exceedance  estimates. 

The  next  section  illustrates  the  computation  of  the  load 
exceedance  functions  with  a flow  chart  and  sample  case. 
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3.  FLOW  CHART  AND  SAMPLE  CALCULATION 


The  calculation  of  the  gust  load  exceedance  parameters  and 
their  use  in  determining  the  design  just  load  are  illustrated  by  a 
flow  chart  and  a numerical  example  in  this  section. 


Output  ; 0 J Input 

parameters  ■ n n>  i parameters 


FLOWCHART  FOR  ESTIMATION  OF  LOADS  ON  WING  AND  VERTICAL 
TAIL  DUE  TO  ATMOSPHERIC  TURBULENCE  USING  A SIMPLIFIED 
POWER  SPECTRAL  TECHNIQUE 
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° VT  CVT 

Calculate 


Vr 


aVTpSVTUK<t>. 


N = — k 
°VT  1TCVT  °VT 

/PN 

FVT=  \t°w  £n 


Sample  Calculation 


Problem:  Estimate  the  loads  on  the  wing  and  vertical  tail  due 

to  turbulence  which  will  not  be  exceeded  more  than 
10  times  in  30,000  hours  of  flight. 
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Necessary  aircraft  parameters  are: 


W = 2950  lb 


S 

w 

cw 

aVT 

SVT 

CVT 


4.88  per  radian 
174.  ft2 

4.57  ft 

it  per  radian 

18.57  ft2 
3.75  ft 


Flight  conditions: 


h = 5000  ft 
U = 219  fps 


The  density  of  air  at  5000  ft  for  standard  temperature  and  pressure 
is 

p = .00205  slug/ft8 
aw  from  Fig.  3 is  6.6  fps 
P from  Fig.  6 is  5 x 10~^ 

For  a rate  of  load  exceedance  less  than  10(5  positive,  5 negative) 
in  30,000  hours,  N the  rate  of  positive  crossings  must  be 

N = ^ = 4.63  x 10"8 

(30,000) (3,600) 
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2W 

U"  ' S,  cw»Sw 


2(295 o: 


23.03 


1500 

337 


1500 

3.75 


Gust  Design  Wing  Load  Factor 


(23,328)  « .43 


kQ  (23,328)  = .066 


a _ 4.88(.00205) (174.) (219.)  ,3  _ Oo78 

AAn 2(295^ 


N = 

o ir  (4 . . 


,066  = 1.007 


. / no->o\  / c.  2N,_  r~5  ( 10  jl.007 

An  = (.  0278)  (6 . 6)  Jin  — g — 

w [_  4 . 62(10”°)  _ 


An  = 1.7  g's 
w ° 

Gust  Design  Vertical  Tail  Load 


K ( oo  400)  - .93 

*VT 

k°\rr  (°°»469)  ~ -017 


by  extrapolation 


= tt(. 00205)18. 57(219.). 93  = 12  2 


- ifr75)<-°w  - °-316 

• 12.2(6.6)in  flli2l_LL  jjlfil 

L 4.62(10'*) 


= 655  lbs 
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